Astronomy & Astrophysics manuscript no. 6240 


© ESO 2008 


February 5, 2008 





Characterization of open cluster remnants * 



O 1 D. B. Pavani 1 ' 2 and E. Bica 2 

O " 



u 1 

Oh: 

<■ 



Instituto de Astronomia, Geofisica e Ciencias Atmosfericas (IAG) Universidade de Sao Paulo - Rua do Matao 1226; 05508-900, Sao Paulo, 
SP; Brazil 

e-mail: daniela@ast.ro . iag . usp . br 

Universidade Federal do Rio Grande do Sul, IF, CP 15051; 912501-970, Porto Alegre, RS; Brazil 
e-mail: bica@if.ufrgs.br 



<-] , ABSTRACT 

' ' Context. Despite progress in the theoretical knowledge of open cluster remnants and the growing search for observational identifications in 
>^ ' recent years, open questions still remain. The methods used to analyze open cluster remnants and criteria to define them as physical systems 



> 

in 



o 



are not homogeneous. In this work we present a systematic method for studying these objects that provides a view of their properties and 
allows their characterization. 

Aims. Eighteen remnant candidates are analyzed by means of photometric and proper motion data. These data provide information on objects 
and their fields. We establish criteria for characterizing open cluster remnants, taking observational uncertainties into account. 
Methods. 2MASS J and H photometry is employed (i) to study structural properties of the objects by means of radial stellar density profiles, 
(ii) to test for any similarity between objects and fields with a statistical comparison method applied to the distributions of stars in the CMDs, 
and (iii) to obtain ages, reddening values, and distances from the CMD, taking an index of isochrone fit into account. The UCAC2 proper 
motions allowed an objective comparison between objects and large solid angle offset fields. 

Results. The objective analysis based on the present methods indicates 13 open-cluster remnants in the sample. Evidence of the presence of 
binary stars is found, as expected for dynamically evolved systems. Finally, we infer possible evolutionary stages among remnants from the 
structure, proper motion, and CMD distributions. The low stellar statistics for individual objects is overcome by means of the construction of 
composite proper motion and CMD diagrams. The distributions of remnants in the composite diagrams resemble the single-star and unresolved 



^— ^ binary star distributions of open clusters. 
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1 . Introduction we concluded that objects older than 1 Gyr amount to w 10 % 

, of the 864 clusters with available ages. 
The study or the space distribution, formation, age, and struc- 

ture of open clusters provides information on their evolution. Qpen duster remnants (here after, OCR) correspond to final 
Open clusters are destroyed over time by the action of both in- evolutionary stages . Numel -ical simulations (Terlevich[l987l de 
ternal forces (mass loss through the dynamical evolution and La Fuente M arcos[l997l[l998l» have shown that the final stellar 
stellar evolution) and external ones, such as interactions with content of QCRs depends Qn ^ masf , function ^ frac . 
the Galactic tidal field^collisions withmolecular clouds, and/or tion of primordial binarieSi Qalactocentric distance, and total 
disc shocking (Friel [T555| Wielen HUB- Portegies Zwart et m The simulations also show that a strong dependence ex- 
al. (EO® confirm central Galaxy tidal-field effects by showing ^ between ±e dissolution timesca i e f an open cluster and 
that star clusters located inside 150 pc of the Galactic center ^ mass evaporation As a conseque nce, OCRs are ex- 
can be dissolved in ~ 50 Myr. In the solar neighborhood, most pected tQ fee rich ^ binary ^ Qwing tQ dynamical procesS es 
open clusters evaporate completely in less than 1 Gyr (Bergond ^ dominate the d i ssolut i on , it is reasonable to assume that the 
et al. [MBonatto et al. EES- However, older objects exist effects of mass segre gation will cause depletion of the low main 
(LyngaE987l Dutra & Bica 2000). By means of the open clus- sequence (MS) in OCRs Considering the above results with 
ter catalogue by Dias et a l. ffl®, updated to January 20041 those of Hsu et al. $199®, Friel ([1993 made further develop- 
ment offprint requests to: D. B. Pavani ments towards an evolutionary scenario. Small and poorly pop- 

* Figures 3,7,9 and, 10 are only available in electronic form via ulated open clusters dissolve primarily from internal dynamical 

http://www.edpsciences.org effects on the order of a few 10 8 Myr, which explains most of 

1 http://www.astro.iag.usp.br/- wilton the open cluster population. Apparently there is no relationship 
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between dissolution rate and location in the Galaxy for low- 
mass open clusters. For richer and more massive open clus- 
ters, the Galactic field effects become significant. Those with 
intermediate masses (~ 500 to 1000M Q ) can survive several 
Gyr, if they are located in the external regions of the disc (e.g., 
NGC 752 and NGC 3680). Open clusters that survive longer 
(e.g. NGC 6791) must be more massive or else located exclu- 
sively in external regions of the disc (NGC 188), or both (Friel 
120021 

From an observational point of view, an OCR can be 
defined as a poorly populated concentration of stars result- 
ing from the dynamical evolution of a more massive system 
(Pavani et al. 2001). Despite theoretical progress, and the grow- 
ing number of observed candidates, open questions still remain, 
such as: (i) Is there a preferential location in the Galaxy for 
OCRs? (ii) Do different evolutionary stages exist among rem- 
nants? (iii) Is it possible to define criteria for the characteriza- 
tion of OCRs? 

Studies in the literature reflect the debate on the obser- 
vational evidence of OCRs. For instance, NGC 6994 (M73) 
was interpreted as a possible remnant of an old open cluster 
by means of high-resolution spectra studies, radial velocities, 
atmospheric parameters, and proper motions (Bassino et al. 
120001 . On the other hand, Carraro (120001 and Odenkirchen & 
Soubiran (2002) interpret NGC 6994 as a random fluctuation of 
the stellar density using CCD photometry in B, V, and I bands 
and high-resolution spectra of the six brightest stars, respec- 
tively. 

Bica et al. ( 120011) present 34 possible open cluster rem- 
nants (POCR) that are located at relatively high galactic lati- 
tudes (\b\ > 15°) and are underpopulated with respect to the 
usual open clusters. They show a significant density contrast 
of brighter stars as compared to the Galactic field. Among the 
objects in that list, some have already been studied in detail 
it was concluded that NGC 1252 is a remnant at d Q = 0.64 
kpc from the Sun, with an age of 3 + 1 Gyr. NGC 1901 is 
a physical system dynamically comparable to the Hyades, at 
d e = 0.45 kpc, and with an age of 0.6 + 0.1 Gyr. Both were 
studied with B,V photometry, from the Tycho-2's proper mo- 
tions and comparisons of object CMDs with model predictions 
for the Galactic field (Pavani et al. l200lT >. Carraro (I20021 con- 
cludes that NGC 7772 and NGC 7036 are remnants. NGC 7772 
is located at d e = 1 kpc, with 14 members and an age of 
1.5 Gyr. NGC 7036 has 17 members and is located at d = 1 
kpc with an age of 3-4 Gyr. He used UBVI photometry for 
color-color diagrams, CMD, and a radial stellar density pro- 
file. NGC 1663 was found to be a remnant of 2 Gyr located 
at d Q - 0.7 kpc (Baume et al. 2003). They employed high- 
resolution spectra of the brightest stars and 2MASS photome- 
try (Skrutskie et al. 1199 7 PI together with Tycho-2's astromet- 
ric data. Pavani et al. (2003) find that Ruprecht3 is a rem- 
nant at d Q = 0.72 kpc with an age of 1.5 + 0.5 Gyr, using J 
and H 2MASS photometry, integrated and individual star spec- 
troscopy, and statistical methods for subtracting field stars from 
the CMD. Bonatto et al. (12004] ) presented NGC 2180 as a prob- 



able missing link between evolved open clusters and remnants 
using J and H 2MASS photometry. They discuss luminosity 
and mass functions, structure, and CMDs. They obtained an 
age of 700 + 70 Myr, d Q = 0.91 + 0.08 kpc, and observed 
mass m b s ~ 47M . Villanova et al. (2004) obtained CCD 
UBVI photometry and medium/high resolution spectroscopy 
for NGC 5385, NGC 2664 and Collinder21 and their fields. 
The analysis using star counts, photometry, radial velocity dis- 
tribution, and Tycho-2's proper motions indicated they are non 
physical objects. Carraro et al. ((2005 ) discuss the nature of 1 1 
possible OCRs by using the Southern Proper Motion (SPM) 
Program 3 combined with 2MASS photometry. They conclude 
that ESO 282 SC26 is a probable physical group. 

The methods employed in the series of studies above are 
not homogeneous. The present study intends to characterize 
OCRs by using a systematic set of observational methods ap- 
plied to a large sample of candidates. We include new POCRs 
and objects in common with previous studies. This is essential 
for establishing criteria for the observational identification and 
determination of OCR properties, taking observational uncer- 
tainties into account. 

In Sect. 2 we present the criteria for object selection and 
their structural properties. In Sect. 3 we describe a statistical 
comparison test between distributions of stars in the CMD. In 
Sect. 4 we discuss the J x (J - H) CMDs and derive ages, red- 
dening values, and distances, together with a classification of 
the objects based on an index of isochrone fitting. In Sect. 5 we 
present the proper motion analysis. In Sect. 6 we characterize 
the OCRs. Finally, concluding remarks are given in Sect. 7. 

2. Methods 

2.1. Object selection 

The selection of POCRs was based on poorly populated stellar 
concentrations included in the catalogues by Alter et al. (1970) 
and Lynga ( 119871 ). The objects have a significant stellar density 
contrast as compared to the background. Among the initially 
selected objects, we kept those presenting evidence of evolu- 
tionary sequences in the CMD. We favored higher Galactic lat- 
itudes to avoid contamination by disc stars. The final sample 
contains 17 objects previously described as star clusters in open 
cluster catalogues and a new one (Object 1), which was found 
on a sky survey plate by one of us (E. B.). Based on DSSH and 
2MASS data we determined accurate equatorial coordinates for 
the POCRs (Tabled. 

Among the 18 sample objects, 8 were previously in- 
cluded in a list of 34 POCRs (Bica et al. |200TT ). Of these, 7 
have been discussed in recent studies: (i) ES0 425SC6 and 
ES0425SC15 were classified as non physical (Carraro et al. 
2005); (ii) NGC 1663 was found to be a remnant (Baume et al. 
120031 ; (iii) Ruprecht 3, NGC 1252 and NGC 1901 were found 
to be remnants (Pavani et al. 120031 Pavani et al. 2001); (iv) 
NGC 6994 (Bassino et al. l2000l Carraro [20001 Ordenkirchen 
& Soubiran 2002 ) is discussed further in the present work. 



2 Two Micron All Sky Survey, All Sky data release, available at 
http://www.ipac.caltech.edu/2mass/releases/allsky 



3 The Digitized Sky Survey available: 

http://cadcwww.dao.nrc.ca/cadcbin/getdss 
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NGC 3680 


OC 


286.76 


16.92 


11:25:38 


-43:14:30 


22.0 ± 1 


NGC 6481 


C 


29.94 


14.93 


17:52:48 


04:10:00 


1.5 ±0.5 


M 73, NGC 6994 


C 


33.71 


-39.93 


20:58:55 


-12:38:03 


4.5 ±0.5 


NGC 6863 


C 


32.27 


-17.99 


20:05:07 


-03:33:20 


1.5 ±0.5 


NGC 1663, OCL-461 


L 


185.88 


-19.66 


04:49:18 


13:09:47 


6.0 ±0.5 


ES0425SC6 


L 


235.39 


-22.28 


06:04:50 


-29:10:59 


3.0 ±0.5 


ES0425SC15 


L 


236.37 


-20.35 


06:14:35 


-29:22:30 


3.0 ±0.5 


Ruprecht 3, OCL-642, ES0426SC33 


C 


238.77 


-14.81 


06:42:07 


-29:27:15 


2.0 ± 0.5 


ES0426SC26 


L 


239.62 


-16.51 


06:36:18 


-30:51:30 


3.0 ± 0.5 


ES0429SC2 


L 


242.59 


-4.16 


07:33:24 


-28:11:15 


5.0 ±0.5 


Ruprecht 3 1 , OCL-7 1 1 , ES068SC9 


C 


250.08 


-5.97 


07:42:58 


-35:35:50 


2.0 ±0.5 


Waterloo 6 


C 


264.82 


-2.67 


08:40:17 


-46:07:55 


1.5 ±0.5 


ESO570SC12 


L 


274.06 


35.89 


11:12:07 


-21:21:00 


5.0 ±0.5 


NGC 1252 


L 


274.08 


-50.83 


03:10:49 


-57:46:00 


6.0 ± 0.5 


NGC 1901, Bok 1, OCL-791.1 


L 


279.03 


-33.60 


05:18:11 


-68:27:00 


5.0 ±0.5 


Object 1 


C 


294.94 


-0.63 


11:41:34 


-62:25:05 


2.0 ±0.5 


ES0132SC14 


C 


308.15 


0.20 


13:36:30 


-62:12:49 


1.5 ±0.5 


ESO383SC10 


L 


312.13 


27.08 


13:31:30 


-35:03:56 


3.0 ±0.5 


Lynga8 , OC1-962, ESQ 226SC1 


C 


333.27 


-0.07 


16:20:04 


-50:13:59 


1.5 ±0.5 



The intermediate-age open cluster NGC 3680 is used as a 
comparison (see also Bica et al. 120011 1. It presents evidence 
of an evolved open cluster dynamical stage and resembles in 
many respects what is expected for OCRs. The central lumi- 
nosity function of NGC3680 is nearly flat from the turnoff to 
the turnover (Bonatto et al. I2004I I. Depletion of the low MS in 
central parts and the presence of a corona rich in low-mass stars 
has been detected in NGC 3680 (see also Anthony-Twarog et 
al. |1997t The 2MASS data are a powerful input in the analyses, 
making possible the study of (i) structural properties, (ii) ob- 
jects and fields by means of a statistical method, (iii) the CMD 
distribution to derive ages, reddening, and distances. 

The kinematical data were extracted from The Second 
U. S. Naval Observatory CCD Astrograph Catalog - UCAC2 
(Zacharias et al. 2004 ). It allows quantitative comparison of the 
proper motions of the object and large solid-angle offset fields. 
The use of complementary methods proved to be essential in 
the present study to deal with completeness effects affecting 
POCRs. 

2.2. Structural properties 

The 18 objects present two characteristic structures in the 
DSS/XDSS images. The first one is poorly populated with 
small angular sizes and is centrally concentrated. Such objects 
are hereafter referred to as compact (C in TableQ}. The second 
structure type is more populated in stars, with larger angular 
sizes and more scattered stars. We refer to them as loose (L in 
TableQ}. 

NGC 3680 is similar to loose POCRs in the present sam- 
ple (Fig. 2, in Bonatto et al. 2004). Figure Q] shows examples 
of R band XDSS images of objects classified by eye as com- 
pact (left panels) and loose ones (right panels). The compact 
objects present a high-density contrast of stars as compared 



to the background and typically include up to 10 stars. Table 
Q] shows the 18 objects of the present sample, together with 
NGC3680. The columns give: (1) designation; (2) structural 
type (C for compact and L for loose); (3) and (4) Galactic lon- 
gitude and latitude, respectively; (5) and (6) J2000.0 equatorial 
coordinates; (7) limiting radius Rn m obtained by means of the 
radial stellar density profile (Figs.|2]and[3j. 

The occurrence of two distinct structures might be associ- 
ated to different evolutionary phases and dynamical processes 
leading to a depletion of stars. Another possibility would be a 
different origin for compact and loose objects related, e.g., to 
formation conditions and position in the Galaxy. 

The 2MASS Point Source Catalogue completeness limi0is 
99% for J < 15.8 and H < 14.1 at \b\ > 30°. Basically, areas 
of incompleteness are at ±75° from the Galactic center, b = 
±1° from the Galactic plane, and radius r < 5° of the Galactic 
center. For these regions the completess limits are J < 14.0 and 
H < 13.5 mag. Figures 7 to 9 of the 2MASS analysis of Release 
Catalog^ allow estimation of completess as a function of 1 and 
b. For the POCRs in Table Q] with \b\ > 14°, the derived limits 
are J < 15.8 and H < 14.1 mag, and for those with \b\ < 6° 
J < 14.0 and H < 13.5 mag. Object 1, ES0132 SC14, and 
Lynga 8 are close to the plane, and the CMD contains very red 
stars, so we applied a color filter (Sects. [4] and [5]) to eliminate 
them, as well as to eliminate the contamination of LMC stars 
in NGC 1901 direction (see Fig [9}. Figures [2] and [3] show the 
radial density profiles for the magnitude-limited extractions for 
the whole sample. 

The radial stellar density profile provides a general view of 
an open cluster structure. We apply this method to POCRs. Star 
counts with 2MASS photometry are carried out in concentric 
rings to take POCR and field areas into account. The extrac- 

4 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_5al.html 

5 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6.html 
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Fig. 1. XDSS R images for 6 sample objects. Left panels: compact objects NGC6481 (top), Object 1 (middle), and Lynga8 
(bottom). Right panels: loose objects (downwards) ES0 425SC15, ES0 426SC26, and ESO570SC12. North is to the top and 
east to the left. 



tions were centered in the equatorial coordinates of Table Q] 
The radial density profiles were obtained from star counts in 
concentric rings with steps in radius depending on the object 
properties, divided by the ring area. The background contribu- 
tion level shown in Figs. [2] and [3] as shaded rectangles corre- 



sponds to the average number of stars present in an outer ring. 
Poisson error bars are significant since POCRs are intrinsically 
poorly populated. 

As comparison, Fig. 4 of Bonatto et al. (120041) shows that 
the radial profile of the evolved open cluster NGC 3680 is de- 
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Fig. 2. Radial stellar density profile of sampled POCR. The average background level is shown (shaded rectangle). Poisson errors 
are also shown. For loose objects the dashed lines show the two-parameter King model fitted to the observed profiles. 



Table 2. Statistical test for POCRs and the open cluster NGC 
3680 



Name 


p 1 

^POCR 


P 




(') 


(%) 


NGC3680 


0.071 


< 1 


NGC 6481 


0.232 


< 1 


NGC 6994 


0.006 


< 22 


NGC 6863 


0.104 


< 11.5 


NGC 1663 


0.017 


< 1 


ES0425 SC6 


0.025 


< 1 


ES0425SC15 


0.027 


< 1 


Ruprecht 3 


0.077 


< 1 


ES0426 SC26 


0.029 


<6 


ES0429 SC2 


0.025 


< 83 


Ruprecht 3 1 


0.197 


< 1.5 


Waterloo 6 


0.184 


< 20 


ESO570 SC12 


0.005 


< 1 


NGC 1252 


0.001 


< 1 


NGC 1901 


0.030 


< 1 


Object 1 


0.459 


< 1.7 


ES0132SC14 


0.380 


< 1 


ESO383SC10 


0.024 


< 1 


Lynga 8 


0.018 


< 1 



scribed by a King model with core radius r c = 2.3 +0.41' and 
a limiting radius R\ lm = 22+1'. Objects in Figs. [2] and [3] il- 
lustrate behaviors in the present sample. Compact POCRs as a 
rule show only one point above the background density. It ba- 
sically corresponds to diameters estimated on DSS and XDSS 
images. Halos are generally absent. King models could not be 



fitted to the profiles of most POCRs. Loose POCRs present 
a density profile with an extended region (possibly a halo) 
characterizing their limiting radius, as illustrated by NGC 1663 
(Fig. The limiting radius arises where an observed profile 
merges with the background. We show the resulting limiting 
radii in Table Q] The loose POCRs NGC 1663, ES0 425SC6, 
ES0 425SC15, ES0 429SC2, ESO383SC10 and NGC 1901 
could be basically fitted by a King model. A comparison with 
the evolved open cluster NGC 3680 (Bonatto et al. |20041 > shows 
that Poisson errors are larger for POCRs. It is important to ex- 
plore loose POCRs as possibly intermediate evolutionary dy- 
namical stages between open clusters and OCRs. Despite the 
loose appearance of ESO570SC12 (Fig. [TJ, its profile resem- 
bles those of compact POCRs, which might also suggest an 
evolutionary transition. 

We conclude that the radial stellar density profiles of the 
sample objects are basically consistent with OCRs that gradu- 
ally lose stars as compared to open clusters, persisting in the 
more internal parts responsible for observed limiting radii. 

3. R 2 statistical test 

Automated procedures for a systematic search of star clusters 
in catalogues have been developed in the past years. Reyle & 
Robin (1200 11 1 combine the density of stars and the integrated 
flux in the K s band. Ivanov et al. (2002) employed the appar- 
ent stellar surface density, K s band luminosity function and the 
distribution of stars in J-K s . Mercer et al. (2005) present an al- 
gorithm that uses statistical methods to locate high stellar den- 
sities aided by criteria of color and magnitude selections. In 
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general the main objective of those searches is to reveal any 
hidden globular clusters, infrared star clusters, and young star 
clusters embedded in HII regions in our Galaxy. For known 
star clusters, there are e.g. studies from Mighell et al. (1996) 
and Kerber & Santiago (2005 ) that applied statistical methods 
to compare populated clusters and fields using the CMD. 

Although OCRs are less populated than open clusters, they 
should still present a stellar density contrast with respect to the 
surrounding field (Sect. l2T2l >. However, it is necessary to adapt 
the open cluster comparison methods owing to the intrinsically 
poorly populated nature of OCRs. Bica et al. (2001) compared 
the density of the bright stars in the POCR area with what is 
predicted for the field from a model of Galactic structure and 
GSQ3 (Lasker et al. 119881 ) counts in an offset equal area. 

We further develop and adapt statistical comparison meth- 
ods to the 2MASS data of POCRs and test the hypothesis 
that the object can be reproduced in radial stellar density and 
CMD distribution by an equal-area field fluctuation. To per- 
form this we employ a nearby large-area offset field, with typ- 
ical angular radius of 30' and separated from the object by 
R\i m + 3 x Rn m . Typically more than 100 sampling field areas of 
size equal to the Ry lm of each object (Table [D were generated. 
Each CMD was divided into boxes with dimensions in color 
A(J — H) = 0.125 and magnitude Aj = 0.5 mag, generating 
a grid. We considered boxes containing at least one star. This 
grid was applied to the object and all samplings in the offset 
field. In each CMD box, the number of stars was counted and 
stellar densities were calculated using the object area. 



6 The HST Guide Star Catalog, Version 1.2, available at http:// 
www-gsss.stsci.edu/ 



Based on the s 2 (Kerber et al. 120011 ), we defined the R 2 
statistics, which reflects the distribution of field fluctuations 
and density contrasts in the CMD plane between field and 
POCRs: 

R 2 = Yj^ -Afield) 2 (1) 

;=i 

where «box is the number of boxes in each grid for POCR (or 
sampled field CMD), p, corresponds to the box density for each 
POCR (or sampled field). Finally, p^aeid corresponds to the box 
density for the total field. 

The ensemble of R 2 for the sampled fields (e.g. Fig. |4]) is 
a measure of the distribution of stars in the CMD of the to- 
tal field, which in turn can be compared to that of the POCR. 
The method proved to be efficient for minimizing undersam- 
pling, which would otherwise be important if we had adopted 
the classical comparison between a populous open cluster and 
equal-area offset field. 

Let jpocr be the position of the POCR R 2 value (Rp 0CE ) 
on the R 2 distribution for iV s sampled fields; we then define the 
probability P of a POCR being a field fluctuation as: 

p _ j _ \jm ~ 7POCrI —\ 

N s /2 

where j m is the median position in that distribution. Table [2] 
gives the P and Rp 0CR values for the POCRs. If the Rp 0CR 
value falls outside the R 2 field distribution, we adopt P < 1%. 
NGC6994, Waterloo 6, and ES0 429SC2 attain P < 22%, 
P < 20%, and P < 83%, respectively, and appear to be field 
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Fig. 4. Examples of R 2 statistical method testing how similar CMDs of objects and respective fields are. Histograms of represen- 
tative field regions within the R 2 interval are given. The arrow shows the object locus. Top panels: NGC 6481 and Object 1 are 
examples of compact objects and appear to be physical ones, while Waterloo 6 appears to be a field fluctuation. Bottom panels: 
NGC 1663, ESO 425 SC6, ESO 426 SC26 are examples of loose objects that appear to be physical ones. 



fluctuations. NGC 6863 presents P < 11.5%, and the remain- 
ing POCRs in Table|2]have very low P values, favoring physical 
systems. 

4. CMD Analyses 

We built 2MASS CMDs of POCRs and offset fields with 
equal extraction radii, corresponding to the region above the 
background field contributions in radial stellar density profile 
(Figs. |2] and 0). We applied a cutoff at J = 15 and excluded 
sources with flags for non-stellar objects for all CMDs in or- 
der to minimize the star field contamination in POCR direc- 
tions. We derived reddening, distance modulus, distance to the 
Sun, Galactocentric distance, and age (Table [3]) with the help 
of solar-metallicity Padova isochrones. In Fig. [5] we present 
the open cluster NGC 3680 core and off-core CMDs. The core 
CMD (left panel) presents a star distribution that is expected for 
OCRs: depletion of the low MS in central parts and the pres- 
ence of a corona rich in low-mass stars. The field of NGC 3680 
at b ~ 17° is uncrowded according to the 2MASS 99% com- 
pleteness limit. NGC 3680 is probably not affected by severe 
crowding, since the brighter limit found in 2MASS crowded 
fieldfO occurs only near the plane or central parts of globular 
clusters. The central MS depletion (Fig. [5]) was also observed 
by Anthony-Twarog et al. (1997). Since POCRs are less pop- 
ulous than open clusters, incompleteness is not expected. The 
brighter stars in the images, which in general characterize a 




J-H 



J-H 



http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6.html 



Fig. 5. Jx(J - H) CMD of NGC 3680 in core (left panel) and off- 
core (right panel) regions. Padova isochrones with solar metal- 
licity and shifted ones according to binarity (dashed line) are 
shown. Photometric errors are indicated. 

POCR (Fig. [TJ, were identified in the CMD, which provided 
constraints for isochrone fittings. 

The NGC 1901 direction is contaminated by LMC stars. We 
eliminated them by means of the color filters shown in Fig [9] 
To verify contamination in the NGC 1901 direction we ana- 
lyzed three CMDs in three filters of fields with the same angu- 
lar size as that of NGC 1901 : (i) the field of the globular cluster 
NGC 1928 projected near the center of the LMC bar; (ii) the 
field of the association of stars NGC 2055 near the east end of 
the LMC bar dominated by a young population, and (iii) the 
field of NGC 2004 with red supergiants. Figures [7] and [9] dot- 
ted lines show the residual LMC contamination in NGC 1901 
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J-H J -H J -H 

Fig. 6. Jx(J - H) CMDs of POCRs. Padova isochrones with solar metallicity, and the same shifted according to binarity (dashed) 
are shown, and photometric errors indicated. 



and offset field directions, respectively. Several POCRs have 
tight MS distributions, such as NGC6481, ES0 425SC6, and 
ES0 426SC26 (Fig.[6]l. A sequence of binaries, together with 
that of single stars in general, encompass most of the POCR 
stars. The sequence of binaries is obtained from that of single 
stars by a shift of Aj = 0.75 mag (Santiago et al. 1996; Elson et 
al. I19981 >. As a rule, offset fields are underpopulated and, scat- 
tered and stars do not show the same brightness range (Figs. 
[8]and|9]l. On the other hand, in the CMDs of Waterloo 6 and 
ESO 132 SC14 (Fig.EJ stars are more scattered. 

We applied a numerical method to quantify the fit of an 
isochrone to a POCR CMD. Stars compatible with single and 
binary isochrones were counted taking error bars into account. 
We obtained their ratio to the total number of stars in the CMD, 
defined as the fitting index //. Best isochrone fits could be 
obtained for all objects, and tentative fundamental parameters 
were derived (Tabled. This suggests that only classical CMD 
analyses may not be enough for conclusive results on the na- 
ture of underpopulated objects. The ensemble of methods will 
settle the more probable physical objects (Sect. [6j. Table [3] 
gives (1) designation, (2) the extraction radius of the data, (3) 
//, (4) age t (Gyr), (5) near IR and (6) optical reddening val- 
ues, (7) distance modulus, (8) distance to the Sun d Q (kpc), (9) 
linear diameter (pc), (10) height from the plane z (kpc), (11) 
Galactocentric distance R gc (kpc), and (12) number of proba- 
ble members. 

Linear diameters were are calculated from angular limiting 
radii (Table [TJ and distances to the Sun. While most POCRs 
have diameters in the range 1-2 pc (TableO, open clusters such 



as M 67 and NGC 3680, at comparable ages and Galactocentric 
distances, have 17.4 and 12.8 pc, respectively (Bonatto et al. 
2004). This suggests that POCRs correspond to open cluster 
cores. The E(B - V) values were coverted from E(J - H) using 
E(B - V) = 3.03 x E(J - H) (Rieke & Lebof sky [19831 1 . The 
number of probable members was estimated from If as was the 
total number of stars in the POCR extraction radius. The num- 
ber of probable members was not estimated for POCRS with 
If < 50%. Typical errors are e(E(J - H)) = 0.02 mag, e(J - 
My) = 0.2 mag, e(d ) = 0.1 kpc, e(z) = 0.015 kpc, e(R gc ) = 
0.4 kpc. The errors combine 2MASS photometric errors with 
CMD fit errors. Using the relation Aj/A v = 0.282 (Rieke & 
Lebof sky [T985I we derived Aj = 2.65 x E(J - H) in the cal- 
culations above. In the sample 13 POCRs with If > 60% are 
described by the respective isochrones well. Three POCRs have 
50% < If < 60%, while Waterloo 6 and ESO 132SC14 have 
// < 50%. 

We determined ages by fitting isochrones to the extracted 
CMDs and using equal area fields as references. No decontam- 
ination procedure was necessary. Most sample POCRs are lo- 
cated at relatively high Galactic latitudes and thus are favored 
by low contamination. The 6 POCRs at \b\ < 6° have small an- 
gular diameters and bluer colors on average than the field stars 
(Fig.0. 

5. Proper motions 

Bica & Bonatto (2005) present internal proper motion distribu- 
tions that are related to binary stars in open clusters. They show 
that proper motions are a useful tool for identifying higher- 
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Fig. 7. Same as Fig.|6]for additional POCRs. 

Table 3. CMDs parameters for POCRs and the open cluster NGC 3680 



Name 


(') 


//(%) 
(%) 


t (Gyr) 
(Gyr) 


E(J - H) 

(mag) 


E(B - V) 
(mag) 


(J - Mj) 
(mag) 


do 
(kpc) 


D 

(pc) 


z 

(kpc) 


R gc 
(kpc) 


pm 


NGC 3680 


3.0 


78 + 18 


1.6 ±0.01" 


0.0 a 


0.0 


10.00" 


1.00 a 


12.8" 


0.30 


8.8 a 


18 


NGC 6481 


1.5 


83 ± 


26 


3.5 


±0.5 


0.06 


0.18 


10.36 


1.18 


1.03 


0.30 


7.04 


10 


NGC 6994 


2.5 


83 ± 


26 


1.0 


±0.2 


0.01 


0.03 


9.77 


0.9 


2.36 


-0.57 


8.0 


10 


NGC 6863 


1.5 


75 ± 


31 


3.5 


±0.5 


0.00 


0.00 


10.40 


1.20 


1.05 


-0.37 


7.06 


7 


NGC 1663 


3.0 


88 ± 


23 


2.8 


±0.2 


0.04 


0.12 


0.20 


1.10 


4.20 


-0.40 


9.13 


15 


ES0425 SC6 


2.5 


71 ± 


18 


2.5 


±0.3 


0.00 


0.00 


10.20 


1.10 


1.92 


-0.41 


8.62 


15 


ES0425 SC15 


3.0 


61 ± 


15 


1.0 


±0.3 


0.01 


0.03 


9.97 


0.99 


1.73 


-0.34 


8.55 


15 


Ruprecht 3 


2.0 


65 ± 


18 


1.5 


±0.5 


0.00 


0.00 


10.20 


1.10 


1.28 


-0.28 


8.60 


13 


ES0426 SC26 


2.5 


70 ± 


17 


1.0 


±0.3 


0.01 


0.04 


10.46 


1.24 


2.16 


-0.35 


8.66 


17 


ES0429 SC2 


1.5 


68 ± 


17 


0.4 


±0.05 


0.03 


0.9 


11.20 


1.67 


4.85 


-0.12 


8.88 


15 


Ruprecht 3 1 


1.5 


78 ± 


24 


1.0 


±0.3 


0.02 


0.06 


9.85 


0.93 


0.81 


-0.10 


8.36 


11 


Waterloo 6 


1.5 


48 ± 


15 


0.2 


±0.1 


0.00 


0.00 


10.5 


1.26 


1.28 


-0.06 


8.21 




ESO570 SC12 


4.0 


62 ± 


20 


0.8 


±0.2 


0.01 


0.03 


9.02 


0.64 


1.86 


0.37 


7.99 


10 


NGC 1252 


4.0 


65 ± 


19 


2.8 


±0.3 


0.00 


0.00 


9.5 


0.79 


2.76 


-0.61 


8.00 


11 


NGC 1901 


5.0 


67 ± 


21 


0.6 


±0.1 


0.01 


0.03 


8.35 


0.46 


1.34 


-0.26 


7.95 


10 


Object 1 


1.5 


50 + 


11 


0.6 


±0.1 


0.00 


0.00 


11.20 


1.74 


2.00 


-0.02 


7.44 


20 


ES0132SC14 


1.5 


33 + 


9 


0.8 


±0.2 


0.00 


0.00 


10.20 


1.10 


0.95 


0.003 


7.37 




ES0 383 SC10 


3.0 


76 ± 


19 


2.0 


±0.5 


0.00 


0.00 


10.00 


1.00 


1.74 


0.45 


7.45 


15 


Lynga 8 


1.5 


56 ± 


19 


2.0 


±0.5 


0.19 


0.57 


10.11 


1.05 


0.92 


-0.001 


7.07 


9 



11 Bonattoetal. (2004). 



velocity stars as unresolved binary cluster members and, as a 
consequence, for mapping and quantifying the binary fraction 
in the CMD. 

Individual POCRs show low stellar statistics as compared 
to open clusters, and composite POCR CMDs might provide 



more insight into their nature. Proper motions extracted from 
UCAC2 allowed a comparison between POCR areas and large 
offset fields. 

We applied a color filter (Bonatto et al. 120051 1 to the CMD 
of a given POCR and offset field stars that removed most back- 
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Fig. 8. Same as Fig.[6]for offset fields. 



ground contamination. The color filter was used to separate the 
CMD of a POCR from background contamination (most disc 
and/or bulge stars), leaving a residual contamination that will 
be taken into account by means of histogram subtractions. In 
the following we employ the modulus of proper motion com- 
ponents 



MP 



aCOsSj 1 +H 2 S . 



(3) 



To be consistent with the CMD analysis, we extracted the data 
for each POCR inside the limiting radius Rn m (Table [T). We 
verified that the correspondence between UCAC2 and 2MASS 
is nearly complete for J < 14.5. 

Histograms of stars in proper motion modulus bins of 2 
mas/yr were built for POCRs with It > 50% and their off- 
set fields. The offset field was scaled to the POCR area. We 
subtracted the offset field histogram from that of the POCR re- 
gion, and the result provides the proper motion distribution in 
Figs. [TOl and fTTl Asymmetries and peaks in the proper motion 
distributions may yield information on the internal kinematics, 
presence of binaries, and possible evolutionary stages. 

In order to analyze intrinsic properties, we transformed 
proper motions into linear velocities on the plane of the sky 
(v p ). The units mas/yr convert to km/s using the distances in 
Tableland 



Vp(km/s) = MPimas/yr) * d(pc) * 4.74 x 10 



(4) 



The objects were divided into two groups (A and B) according 
to their v p distributions characterized by peaks or flatter distri- 
butions. The observed median velocities (vm) were derived for 



rest velocity corrections (Table @), and the median is suitable 
for non Gaussian distributions. 

Finally, composite velocity histograms were built: (i) vm in 
each POCR histogram was used to correct the diagram to rest 
velocity, providing it in the form (v p — Vm)', (ii) Group A is char- 
acterized by a well-defined low-velocity peak and is formed by 
the loose POCRs NGC 1663, ES0425 SC15, ES0426 SC26, 
NGC 1 252, NGC 1 90 1 , and the compact one Ruprecht 3 1 . They 
have (vm < 71 km/s). Group B is characterized by higher- 
velocities, 75 < vm < 140 km/s, and is formed by the 
compact objects NGC 6481, NGC 6863, Ruprecht 3, Object 1, 
Lynga8 and the loose ones, ES0 425 SC6, ESO570SC12, and 
ESO383SC10. 

Figure [12] presents the composite histograms (v p - vm) for 
Groups A and B. Group A shows a well-defined low-velocity 
peak. Group B shows a double peak. Table |4]gives the POCRs 
belonging to each group. It is outstanding that the (v p - vm) 
composite histogram reveals properties of POCRs only seen 
before in open clusters (Bica & Bonatto 2005), which can- 
not be observed in individual POCRs owing to low statistics. 
The double-peak distribution in Group B resembles that of sin- 
gle (low velocity) and unresolved binary stars (high velocity) 
as observed in off-core regions of M67 (e.g. Fig. 5 in Bica & 
Bonatto 120051 

Several objects were not included in the group analysis be- 
cause of doubt about their nature or they are field fluctuations 
because the indices are indicating non physical objects (Table 
|5). Although ES0 429SC2 has a proper motion distribution 
compatible with Group A and has If = 68 + 16%, it was not in- 
cluded in the composite histogram due to the probability of be- 
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Fig. 9. Same as Fig.[6]for additional offset fields. The thick line in the NGC 1901 CMD indicates contamination of LMC stars 
eliminated by color filters. 



Table 4. The groups and median velocity values 



Name 


Type 


v M (km/s) 


Group A 


NGC 1663 


L 


68 


ES0425SC15 


L 


47 


ES0426 SC26 


L 


70 


Ruprecht 3 1 


C 


44 


NGC 1252 


L 


68 


NGC 1901 


L 


26 


Group B 


NGC 6481 


C 


78 


NGC 6863 


C 


137 


ES0425 SC6 


L 


83 


Ruprecht 3 


C 


83 


ESO570 SC12 


L 


109 


Object 1 


C 


82 


ESO383SC10 


L 


113 


Lynga 8 


C 


75 



ing a field fluctuation (Table 0. The CMD of NGC 6994 (Fig. 
[6]l shows a gap of «3 mag. between bright and faint stars that 
is improbable in a turnoff (TO) to low MS distribution. The 
statistical method (Table [U also favors a non physical nature. 
These results point to NGC 6994 being a field fluctuation, in 
agreement with Carrara (2005 ) and Odenkirchen & Soubiran 
(120021 . 

Composite open cluster CMDs have been a tool employed 
since Sandage ( 119581 and were useful for modeling stellar evo- 



lution. In the present study we applied this principle to OCRs 
and POCRs. By means of reddening and distance corrections, 
we obtained absolute CMDs, which were in turn combined to 
build composite CMDs for Groups A and B, improving the stel- 
lar statistics in the diagrams (Fig.[T3~ll. Considering the (v p -vm) 
distributions for Groups A and B, stars with velocities higher 
than 25 km/s are more probably related to unresolved binary 
stars, based on analysis of the populous open cluster M 67 (Bica 
& Bonatto 2005). We plotted those stars in the composite CMD 
and they tend to be located to the right of the MS of single stars, 
as expected (Fig.Qj). 

Group A (Fig.Qj) shows the low-velocity stars distributed 
along the MS for different ages, while the higher-velocity stars 
are biased towards the right side of the isochrones. This is the 
locus encompassed by the color and magnitude limits expected 
for binary stars. In the composite CMD of Group B (right 
panel) a similar behavior is seen, and stars around the 3.5 Gyr 
isochrone suggest that higher-velocity stars tend to replace the 
single-star MS. This loss of single-stars to the field would tend 
to cause a higher fraction of binaries in such OCRs, in agree- 
ment with the results of N-body simulations (de La Fuente 
Marcos [l997l [19981 . 

The composite diagrams of Groups A and B suggest pos- 
sible evolutionary stages among the studied remnants, where 
Group B would be more dynamically evolved than Group A, 
and in general loose objects would be less evolved than com- 
pact ones. 



12 



N 



N 



N 



2 
10 
8 
6 
4 
2 



A' 



n i 1 r 

NGC 1663 



m 



+ 




ESQ 425 SC15. 



7 



+ 



+ 



NGC 1901 



i 1 i 1 r 



\ — — V 



iH 



h n 



10 



20 30 
PM (mas/yr) 



40 



10 



20 30 
PM (mas/yr) 



40 



50 



Fig. 10. Left panel: we illustrate the modulus of the proper motion distribution for the object region (histogram) and offset field 
(shaded area). Right panel: modulus of the proper motion distribution after field subtraction: these POCRs are characterized by a 
well-defined peak. 



6. Discussion 

The present set of methods for the study of POCRs provided 
a view of their properties and criteria for their characteriza- 
tion. Of the 18 sample objects, ages, reddening values, and dis- 
tances could be determined for 16 of them having If > 50% 
(Table |3). In Table [5] we present a final POCR classification 
according to different criteria: (i) radial stellar density profile 
considering Poisson errors, (ii) probability of less than 15% for 
the POCR CMD to be representative of the offset field CMD, 
(iii) homogeneous distribution of stars in the CMD and fitting 
index If > 55%. Columns give the partial classifications in 
each method, where "y" stands for the attained minimum value, 
while "n" stands for those not attained. The final results are in- 
dicated in the last column, where "OCR" is an object we con- 
cluded is an open cluster remnant, FF a field fluctuation, and 
POCR is maintained as a possible open cluster remnant. 

Structural and CMD analyses together are powerful tools 
for inferring the nature of POCRs. As an example, although 
NGC 6994 has a reasonable density profile, it has heteroge- 
neous distribution of stars in the CMD (Fig. |6j, strongly sug- 
gesting a non-physical object. 

The statistical method comparing object and field CMDs 
has provided complementary results. It tests whether an object 



has a CMD compatible with that of the field. We conclude that 
72 % of the studied objects have a probability p < 2% of being 
representative of the field (Table |2J. 

Proper motions are indicators of the internal kinematics of 
the studied objects. Since we are dealing with intrinsically un- 
derpopulated stellar systems, it is important to build compos- 
ite CMDs and proper-motion diagrams. The OCR composite 
CMDs show similarities with those of individual open clus- 
ters. The OCR composite proper-motion diagrams suggest their 
classification into two distinct patterns: Group A with one and 
Group B with two pronounced velocity peaks. The distinction 
between binary and single stars and their mapping in the CMD 
suggest evolutionary stages. 

Figure [14] shows the age distribution of OCRs and POCRs 
with Ij > 50% separated in Groups A and B. The peak of the 
distribution is in the range 0.5-1.0 Gyr, but older objects also 
exist. Figure Q3] shows the objects in Groups A and B in terms 
of age as a function of Galactocentric distance R gc (left panel) 
and height from the plane \z\ as a function of R gc (right panel). 
For comparison we included 230 open clusters from Chen et 
al. (120031 . The Sun is assumed to be at 8 kpc (Reid [19931 . 
Despite the small sample and the error bars, if we consider that 
the Group B objects tend to be located inside the Solar cir- 
cle, this configuration could suggest that tidal effects related 
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Fig. 11. Same as Fig[10l for POCRs with a flatter distribution of the modulus of proper motion. 



to R gc affect the evolution of remnants (Table [3]). Considering 
a shorter distance (R gc = 7.2 kpc) of the Sun to the Galactic 
center (Bica et al. 120061 . tidal effects would be stronger. 

7. Concluding remarks 

We studied 18 POCRs by means of 2MASS J and H pho- 
tometry and UCAC2 proper motions. We analyzed structural 
properties and distributions of stars in CMDs of POCRs and 
offset fields using a statistical method. Ages were derived, to- 
gether with reddening values and distances to the Sun, using 
the CMDs and Padova isochrones. We defined an index of 
isochrone fit. All these methods are useful criteria for remnant 
characterization. The analysis indicates the presence of 13 open 
cluster remnants in the sample (Table |5). 

Open clusters are known to contain a significant fraction 
of multiple systems as a consequence of their dynamical evo- 
lution. The fraction of multiple systems increases and tends to 
concentrate in the central regions, thus changing the initial spa- 
tial distribution of stars (Takahashi & Portegies Z wart [20001. 
Numerical simulations show that open cluster remnants are ex- 
pected to be rich in binary stars. As an example, de La Fuente 
Marcos (1997) stopped simulations when the cluster popula- 
tion was 10 stars and the number of the remaining binaries 
was 3 or 4. Recent studies comparing observational and the- 
oretical results indicate that cluster remnants can form highly 
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Fig. 12. Distribution of velocities projected on the sky reduced 
to the rest velocity. Left panel: Group A. Right panel: Group B. 
A velocity bin of 20 km/s was used. 

hierarchical high-order multiple systems; and considering that 
the number of such objects can be larger, cluster remnants may 
be sufficient to account for all observed higher-order multiple 
systems (Goodwin & Kroupa 2005). 

The present overall results indicate evolutionary differences 
among sample objects. POCRs show two distinct structures, 
compact and loose, and we verify a correspondence between 
the radial stellar density profiles and these structures. Sample 
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Table 5. Final classification for the POCRS 



Nome 


Type 


Profile 


r" 


j 


Q 


NGC 3680 


L 


v 

y 


v 

y 


v 

y 


oc 


NGC 6481 


C 


y 


y 


y 


OCR 


NGC 6994 


C 


y 


n 


n 


FF 


NGC 6863 


C 


y 


y 


y 


OCR 


NGC 1663 


L 


y 


y 


y 


OCR 


ES0425 SC6 


L 


y 


y 


y 


OCR 


ES0425 SC 15 


L 


y 


y 


y 


OCR 


Ruprecht 3 


C 


y 


y 


y 


OCR 


DOU 4ZO JL ZD 


T 


y 


y 


y 


ULK 


ES0429SC2 


L 


y 


n 


y 


POCR 


Ruprecht 3 1 


C 


y 


y 


y 


OCR 


Waterloo 6 


C 


y 


n 


n 


FF 


ESO570SC 12 


L 


y 


y 


y 


OCR 


NGC 1252 


L 


y 


y 


y 


OCR 


NGC 1901 


L 


y 


y 


y 


OCR 


Object 1 


C 


y 


y 


n 


POCR 
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Fig. 13. Absolute CMDs for Groups A and B (left and right panels, respectively). Stars with velocities (v p - v«) higher than 25 
km/s (filled circles) are superimposed on the CMDs. Padova isochrones are shown as reference. 



objects are basically consistent with OCRs that gradually lose 
stars as compared to open clusters. The composite proper mo- 
tion diagrams shows evidence of binary stars, as expected from 
high binary fractions in dynamically evolved systems. The 
presence of a sequence of binaries in the composite CMD, re- 
placing the MS of single-stars in Group B, suggests the exis- 
tence of distinct stages along the dynamical evolution of open 
cluster remnants. In addition, most of the compact objects are 
included in Group B and tend to be more internally located in 
the disc. If two distinct structures might be associated to dif- 
ferent evolutionary phases, then transitional objects should be 
detected. ESO 570 SC12 is visually loose (Fig.[TJ but was com- 



pact in the radial density profile (Fig. [3). However, a large sam- 
ple would be necessary to constrain such possible links. For the 
moment we cannot exclude the possibility that Groups A and 
B might be descendants of different types of open clusters. 

For the first time a study presents a considerable number 
of POCRs analyzed by means of the same methods for the 
characterization of their common properties. We verify that is 
necessary to combine different methods in order to constrain 
the object nature owing to data incompleteness, observational 
uncertainties, and intrinsically low statistics. The existence of 
compact and loose POCRs suggests relatively stable remnant 
structures that dynamically survive beyond open clusters. 
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Fig. 14. Age histogram for OCRs and POCRs. 

We confirm NGC 1252, NGC 1901, Ruprecht3, and 
NGC 1663 as open cluster remnants, and NGC 6994 as a 
field fluctuation, in agreement with previous studies (Pavani 
et al. I200T1 Pavani et al. [20031 Baume et al. 120031 Carraro 
2000 1 Odenkirchen & Soubiran 120021) . Evidence is found 
that NGC 6481, NGC 6863, ES0425SC6, ES0 425SC15, 
ES0 426SC26, ESO570SC12, Ruprecht31, ESO383SC10, 
and Lynga 8 are also open cluster remnants, while Waterloo 6 is 
a field fluctuation. ES0429 SC2, Object 1, and ESO 132SC14 
are maintained as POCRs for more detailed future studies. 

The present work offers observational evidence of the ex- 
istence in nature of open cluster remnants, as expected from 
theoretical and numerical studies. As prospective work to an- 
swer fundamental questions, we mention (i) would the distinc- 
tion between compact and loose structures be related to that of 
the progenitor open cluster and/or the evolutionary processes 
and (ii), if evolutionary stages in remnants are confirmed, what 
detailed processes would cause them? 
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